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STANDARD DEVIATION OF VERTICAL WIND SHEAR IN THE LOWER ATMOSPHERE 


SUMMARY 


The standard deviation a of vertical two-point longitudinal veloc- 
ity fluctuation differences is analyzed experimentally with eleven sets 
of turbulence measurements obtained at the NASA 150-meter ground winds 
tower site at Cape Kennedy, Florida. It is concluded that the a! u* Q is 
proportional to (fAz/u^o) 1 ' 3 , where the coefficient of proportionality 
is a function of f z/u* 0 and u*o/fL 0 . The quantities f and L 0 denote 
the Coriolis parameter and the surface Monin-Obukhov stability length, 
respectively; u* 0 is the surface friction velocity; Az is the vertical 
distance between the two points over which the velocity difference is 
calculated; and z is the mean height of the mid-point of the interval 
Az above natural grade. The results of the analysis are valid for 
20 < -u* 0 /fL 0 < 2000. 


I. INTRODUCTION 


Statistical estimates of wind shear in the planetary boundary layer 
are important in the design of V/STOL aircraft (vertical/short take-off 
and landing aircraft). The STOL aircraft are sensitive to variations 
in lift forces resulting from vertical variations of the horizontal 
wind during take-off and particularly during landing operations. VTOL 
aircraft control systems must be designed to counteract extreme tipping 
moments resulting from vertical variations of the horizontal wind for 
specified acceptable values of risk of exceeding the design wind shear 
conditions. Similarly, wind shear information will also be necessary 
in the design of the forthcoming NASA Space Shuttle for both the launch 
and landing phases. During the launch phase the Shuttle will experience 
tipping moments; during the landing phase, it will experience fluctuating 
lift forces. The vertical two-point longitudinal velocity differences 
that will be experienced by these vehicles can be divided into two parts: 
the wind shear resulting from the mean flow, and the fluctuating or random 
part resulting from turbulence. The mean shear can be calculated with 
currently available models of the mean flow wind profile in the atmos- 
pheric boundary layer [1,2, and 3], However, the fluctuating part must 
be specified statistically for engineering problems, and this requires 
knowledge about the statistical moments or the distribution function of 
the shear. A great deal of information can be gained by studying the 
second moment of wind shear. In addition, the square root of the second 



moment, the standard deviation, of a random variable is used as a 
standardizing parameter to render the random variable nondimens ional. 
Accordingly, the analysis in this report is restricted to an examina- 
tion of the second moment of vertical two-point longitudinal velocity 
differences, or, in other words, wind shear. 

It is hypothesized that the standard deviation of the wind shear 
scaled with the surface friction velocity u* D is a universal function of 
the nondimens ional parameters fAz/u* 0 , fz/u* 0 , and u* 0 /fl 0 > where f is 
the Coriolis parameter, L Q is the surface Monin-Obukhov length, Az is 
the vertical distance between the two points over which the shear is cal- 
culated, and z is the mean height of the mid-point of the shear interval 
Az above natural grade. 

From the eleven cases of turbulence in unstable air selected for 
analysis, it is concluded that cr/u* 0 is indeed a function of the non- 
dimensional parameters listed above. Empirical mathematical expressions 
are obtained to represent this function. It appears that a/ u* 0 is pro- 
portional to (fAz/u-fro) 1 ' 3 , the coefficient of proportionality being a 
function of fz/u, vo and u* 0 /fL 0> The results of the analysis are valid 
for 20 < -u* 0 /fL 0 < 2000. 

The author expresses his thanks to Mr. Archie Jackson and Mrs . Ella 
Mae McAllister of the Marshall Space Flight Center's Computation Labora- 
tory. Mr. Jackson developed the program for the calculation of the two- 
point standard deviations. Mrs. McAllister developed the remaining 
programs which were used to calculate the non-dimensional quantities. 

In addition, thanks go to Mr. Julian Nelson and Mr. Douglas Mackiernan 
of the Aero-Astrodynamics Laboratory for their help in the preparation 
of this report and the analysis of the data. 


II. THE DATA SOURCE 


The data analyzed in this report consist of eleven sets of longi- 
tudinal turbulent velocity fluctuation time histories digitized at 
0.2-second intervals with 18,000 data points per time history. The 
longitudinal velocity fluctuations were calculated with horizontal wind 
speed and direction data measured at the 18- , 30- , 60- , 90- , 120- , and 
150-meter levels at the NASA meteorological tower site at the Kennedy 
Space Center with Climet (model Cl-14) wind sensors. The measurements 
were taken during the daytime in unstable air. The computation procedures 
used to calculate the longitudinal turbulent velocity fluctuations are 
discussed by Fichtl and McVehil in reference 4. Supporting temperature 
measurements were made at the 18- and 30-meter levels with Climet (model 
Cl-016) aspirated thermocouples. (Details concerning the instrumentation 
at the NASA tower site can be found in a report by Kaufman and Keene [5]; 
the surface roughness lengths associated with this site are discussed by 
Fichtl and McVehil in reference 4.) 
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III. THEORETICAL CONSIDERATIONS 


The instantaneous longitudinal velocity u(x,y,z,t) in a horizontally 
homogeneous atmospheric boundary layer can be represented as 


u(x,y,z,t) = u(z) + u* (x,y,z,t), 


( 1 ) 


where u(z) is the temporal mean wind at height z, and u ! (x,y,z,t) is a 
longitudinal turbulent velocity fluctuation. Upon evaluating equation 
(1) at heights z x and z 2 (z 2 > z x ) for fixed values of x,y, and t, we 
may calculate the two-point velocity difference 


Au(z 1 ,z s ,t) - ^(z^za) + Au' (zi,z 2 ,t), (2) 

where 

= u(z 2 ,t) - u(z x ,t) (3) 

Au(z x ,z 2 ) = u(z 2 ) - u(z x ) (4) 

£u' (z x ,z 2 ,t) = u' (z 2 ,t) - u' (z x ,t). (5) 


The dependence of Au and Au* on x and y is understood. For engineering 
purposes, the quantity Au(z l9 z 2 ) can be treated as a known quantity 
that can be calculated with a power law wind profile like 

u(z) = u(z 1 )(z/z 1 ) a , (6) 


where a is a positive constant, or with the more elegant wind profile 
formulations for the barotropic and baroclinic Ekman planetary boundary 
layers by Blackadar [1], For values of z 2 < 30 m, the mean wind profile 
law appropriate for a Monin layer can be used: 


u(z) = (in - i|r(z/L 0 )^, (7) 

where u* 0 and L Q denote the surface friction velocity and Monin-Obukhov 
stability length, z Q is the surface roughness length, \a ± is von Harman' s 
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constant with numerical value approximately equal to 0.4 and \|r(z/L 0 ) is 
a universal function of z/L q [3], On the other hand, the quantity 
Z\u f (z,z 2 ,t) is a rapidly varying random function of time, and thus must 
be treated statistically. 

The longitudinal time histories were differenced in the vertical to 
yield fifteen longitudinal velocity difference time histories for each 
case. Unbiased estimates of the second moment were calculated with the 
formula 


n 


i=l 


( 8 ) 


where a is the standard deviation of the random longitudinal shear and 
n is the number of data points used in the calculation. 

IV. DIMENSIONAL ANALYSIS 


The Monin-Obukhov similarity hypothesis [3] predicts that the 
standard deviation of two-point longitudinal velocity differences, scaled 
with the surface friction velocity u, vo , are universal functions of z^/Lq 
and z 2 /L d or z/L 0 and Az/L d , where 


z 


Z X + z 2 
2 


Az = z 2 - Zj 


(9) 


Thus, the Monin-Obukhov similarity hypothesis leads to the conclusion 
that 


o ! u * 0 = F(z/L q , Az/L d ), 


( 10 ) 


where F is a universal function of z/L 0 and Az/L q . However, as one pro- 
ceeds upward out of the Monin layer into the Ekman layer, it is reason- 
able to hypothesize the cr/u*o is a function of z/L Q , Az/L 0 and other 
parameters that characterize the action of Coriolis forces and baroclinic 
effects. Blackadar and Tennekes [6] have shown that, in the horizontally 
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homogeneous, neutral, barotropic Ekman layer, u^/f is the appropriate 
length scale, where f is the Coriolis parameter. This result was pro- 
duced by analyzing the mean flow momentum conservation equations along 
with the turbulent energy equation. Thus, as a working hypothesis, one 
might assume that 0 / u^ Q is a function of the independent variables 


£ 2 , z, L 0 > u* Q /f. 


(U) 


According to Buckingham 1 s theorem [7], the number of independent non- 
dimensional quantities that can be constructed from this list of 
variables is three. Now there are eighteen possible ways to select 
this group of three independent nondimens ional quantities such that z 
and Az appear within each group. However, any one of the possible 
groups can be derived from the remaining seventeen groups. Accordingly, 
we are free to select any one group we like to perform a data analysis; 
thus, we hypothesize that 




£ z U * qN \ 

u *o’ f V 


( 12 ) 


This hypothesis includes the dominant effects in the Ekman layer: 
atmospheric stability, and the action of Coriolis forces in the hori- 
zontal mean flow momentum conservation equations. It does not include 
baroclinic effects resulting from a height-dependent horizontal pressure- 
gradient force. Nevertheless, the hypothesis given by equation (12) 
seems to be appropriate in view of the fact that the outstanding feature 
of the Ekman layer is the action of Coriolis forces and that inclusion 
of baroclinic effects will only lead to refinements. This hypothesis 
will be used to analyze the standard deviation of vertical two-point 
longitudinal velocity fluctuation differences. The procedures used in 
the calculation of u* 0 and L Q are given in appendix A. The data used in 
these calculations and the results of the calculations are given in 
appendix B. 


V. DATA ANALYSIS 


To estimate the universal function F on the right-hand side of 
equation (12), the nondimens ional standard deviation a/u* 0 was plotted 
as a function of fAz/u* 0 for various categories of fz/u* D and u* 0 /fL 0 . 
Each case is characterized by a single value of u* 0 /fl j 0 ; however, within 
each case, the quantity fz/u* 0 can vary. 
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Figure 1 is a plot of a/u, vo as a function of f Az/u* 0 for the 
various categories of fz/u AO and u* 0 /fL 0 . The data appear to show that, 
for each of the selected categories, a/u* 0 is proportional to (fAz/u^) 1 / 3 . 
In addition, within each category of fz/u* 0 , the values of a/u* 0 tend to 
sort out according to the parameter u* 0 /fL 0 ; i.e., for fixed values fz/u^ 0 
and fAz/Uko> the quantity a/u*. 0 tends to be a decreasing function of 
-u* 0 /fL 0 . There also appears to be a dependence of a/u* Q on fz/u^o; i.e., 
for fixed values of fAz/u* 0 and u* 0 /fL 0 , the quantity a /u i(Q is a decreas- 
ing function fz/u* Q . 

To see these effects better, it was assumed that a /u^ 0 is proportional 
to (fAz/Ufco) 1 / 3 for all categories of the remaining independent variables, 
so that 


fz u *o\ /fAz/ u* D v l/3 

, u *o ’ f V V 0 - 001 ) ’ 


(13) 


where G is a universal function of ,fz/u^ 0 and u* 0 /fL 0 . The function G 
was determined by fitting equation (13) to the data in figure 1 for each 
category of fz/u* 0 and u*o/fL>o. The results of these calculations are 
shown in figure 2. We see from this figure that G and thus cr/u^ 0 are 
clearly decreasing functions of fz/u^ Q and increasing functions of u* 0 /fL 0 . 

To give the experimental results mathematical expression, we assume 
that the curves in figure 2 are shape-invariant with respect to variation 
in u* 0 /fL 0 and that translation only with respect to the G-axis is required 
to produce a collapse of these curves. This means that there exists a 
function H(u* 0 / fL o) suc h that 


G 





(14) 


In view of the experimental results in figure 2, these assumptions appear 
to be reasonable. To calculate the function H(u* 0 /fL 0 ), we subtract the 
experimental estimate of G(fz/u* 0 > -150) from each curve in figure 2 and 
then average each of the resulting difference functions over the 
experimental range of variation fz/u* 0 . The experimental esti- 
mates of H(u* 0 /fL 0 ) for the various values of u* 0 /fL 0 are given in 
figure 3. The function 


H( u * 0 / fL 0 ) - -0.072 In 


u *o/ fL o\ 
150 ) 


(15) 
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appears to summarize these results reasonably well. To improve our esti- 
mate of G(fz/u* 0 , -150), we subtract H(u^ 0 /fLg) as given by equation (15) 
from the experimental values of G(fz/u* 0 , u* 0 /fL 0 ) given in figure 2 and 
average the resulting difference functions over the experimental range 
of variation of u* 0 /fL Q for each category of fz/u* 0 . The revised esti- 
mate of the function G(fi/u* 0 , -150) is depicted in figure 4, and the 
analytical expression 


G ( f z/u* 0 » 


-150) = 1 + 2.77 (psr) 


(16) 


summarizes these data. The results in figure 4 show that G(fz/u^ 0 , 
-150) -> 1 as fz/u* 0 -» This means that a/u* Q becomes asymptotically 
height- invariant as fz/u* 0 oo. The function (16) contains this 
asymptotic behavior. 

1. 

The combination of equations (13) through (15) yields 




0.072 In 


150 /J \ 0.001 


(17) 


According to equation (17) a/u* 0 becomes large as fz/u* G 0 and 
fAz/u* 0 oo. This means that a becomes large as the height i decreases 
and the interval 4sZ increases. In addition, a/ u* D decreases as -u* Q / fL Q 
increases. However, equation (17) predicts that a/u*o < 0 at sufficiently 
large values of -u* 0 /fL 0 . This is nonsense and is an apparent shortcoming 
of the empirical formula (17). Negative values of a/u* D are predicted 
with (17) when 


fz/u* 0 > 0.001 


L 

0.072 In 


2.77 


u *o /fL o 

150 



(18) 


provided -u* D / fL Q > 4 x 10 s . However, values of -u* 0 / FL 0 > 4 x 10 s are 
very rare if we exclude the region of the earth on and about the equator. 
A physically realizable upper bound value on -u*o/fL 0 in unstable air is 
probably on the order of 10,000 if we exclude the equatorial region of 
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the earth. Thus, the shortcoming of predicting negative values of 0 / u* Q 
with equation (17) is only apparent and would not arise in practical 
situations. The function H(u^ 0 /fL 0 ) probably tends to approach a con- 
stant value as -u* 0 /fL Q becomes large and if this occurs in the interval 
1,000 < -u* 0 /fL 0 < 10,000, then the error in the estimate of o/u* 0 
resulting from (17) is approximately 5 percent. As -u* 0 / £L 0 -» 0 from 
positive values of -u* 0 /fL 0 , a/ u, vo oo. The function H(u* 0 /fL 0 ) should 
approach a finite value as -u* 0 /fL 0 0 because of the finite nature of 

a. Accordingly, H(u* 0 /ff>o) does not have logarithmic behavior as 
-Ufc 0 /fL 0 -» 0 as predicted by (15). The precise value of -u* 0 /fL 0 where 
H(u* 0 /fL 0 ) departs from the behavior stated by (15) can only be deter- 
mined with measurements of 0 / u, vo associated with values of -u* 0 /fL 0 < 20. 
In view of these comments, it is recommended that (17) be used to obtain 
estimates of a only for those situations in which u* 0 / fL 0 is contained 
within the experimental range of variation in this report, i.e., 

20 S -u^ 0 /fL 0 § 2,000, Estimates of a/u* Q can be obtained with equa- 
tion (17) for 2,000 < -u* 0 /fL 0 < 10,000 which probably have errors equal 
to 5 percent. 

The correlation coefficient of the longitudinal wind fluctuations 
at z x and z 2 is 


R(u' (z^.u' (z 2 )) 


^(Zi) + c£( Z 2 ) - a 2 (z 1} z 2 ) 
2a u (zi) a u (z 2 ) 


(19) 


where a u (z) is the standard deviation of the longitudinal wind fluctua- 
tions at height z. Fichtl and McVehil [4] find that for the unstable 
boundary layer 


— « 1.9(z/18)'°- 14 . (20) 

u *o 


This result was obtained by effectively averaging together measurements 
of cr u /u* 0 obtained during various degrees of instability in the first 
150 m of the atmospheric boundary layer at Cape Kennedy, Florida. How- 
ever, the recent results of Panofsky and Mirabella [9] appear to show 
that cr u /u* 0 is relatively insensitive to changes in stability and thus 
u* 0 /fL 0 . Now, (20) shows that a u / u, vo is a slowly varying function of z 
so that we will assume a u (z x ) ^ 0 u (z 2 ) and set 0 u /u, vo = 1.5. This value 
of 0 / UfcQ occurs at z ~ 80 m. Thus, we may write (19) as 

( 0 /u* o ) 2 

R(u' (zi),u' (z 2 )) - 1 - — ^5 — . (21) 
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We may conclude from (17) and (21) that in the unstable boundary layer 
the interlevel correlation coefficient of the longitudinal velocity 
fluctuations increases as the instability of the boundary layer and 
fz/u* 0 increase and as fAz/u* 0 decreases. 


VI. CONCLUDING COMMENTS 


1. We have experimentally analyzed the standard deviation of 
vertical two-point longitudinal velocity fluctuation differences, and 
have found that cr/u* 0 is indeed a function of fz/u^Q, fAz/u*p> and 
u*o/fL 0 . In particular, cr/u* 0 is proportional to (fAz/u^) 1 / 3 . 

Although the analysis did not include baroclinic effects, we believe 
them to be of second order. The fact that a 2 is proportional to 
(Az) 2 / 3 is, in all likelihood, related to the fact that the structure 
function of longitudinal velocity fluctuations at a given height is 
proportional to (Z&) 2 / 3 for sufficiently small values of the horizontal 
space lag Ax when an inertial subrange is present. The precise nature 
of this relationship is not obvious; however, it will be the subject 

of future investigation by the author. 

2, This study is the first part of a four-part study to define 
the statistical properties of fluctuating wind shears in the atmospheric 
boundary layer. The second part of the study will be concerned with 
determining the dependence of the third and fourth moments, p, 3 and [i 4 , 
on fAz/u* 0 , fz/u^ 0 , and u^ 0 /fL 0 in unstable air. The third part will 

be concerned with determining a class of distribution functions which con- 
tain four parameters that will adequately describe the observed wind 
shear distribution functions in unstable air. Preliminary results 
appear to show that the Pearson system of distribution functions can 
be used; to represent the data. The four parameters which will appear 
in the distribution function will be completely specified upon specify- 
ing the s first four moments of the wind shear. The data in this report 
will be used in parts two and three of the total study. The fourth 
part will be concerned with extending the results of parts one, two, 
and three to (1) the neutral boundary layer with turbulence data for 
relatively high wind speed conditions (u(18m) § 15 m sec -1 , say) and 

(2) stable boundary layers associated with nighttime conditions. 
Accordingly, the overall goal is to develop models of a, p, 3 and p, 4 
such that these parameters can be calculated for given values of the 
height parameters z and Az once the external parameters L 0 , u* Q , and f 
are specified. Substitution of a, p 3 , |_i 4 and the mean wind shear (see 
section III) into the distribution function of wind shear will facili- 
tate the calculation of wind shear for any percentile level of occur- 
rence for a given state of the boundary layer. Furthermore, the model 
will permit climatological studies of wind shear with rather straight- 
forward statistical procedures for specified locations upon specifying 
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the climatological statistics of L 0 and u* D (joint distribution function). 
The distribution function of these external parameters could be extracted 
from standard meteorological observations. 
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Figure 1, Experimental values of a/u* Q for various categories 

of fz/u* 0 and -u* 0 /fL_ 


O a < OO 




Figure 2. Estimate of the function G(f5/u. 
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Estimate of the function H(u* 0 /fL 0 ) 






APPENDIX A 


Calculations of the Scaling Parameters and L Q 


In the unstable Monin layer, the dimensionless mean flow 3hear is a 
universal function of z/L 0 , so that 



i- V'V. 


(A-l) 


where u(z) is the mean wind speed at height z, is von Karman' s constant 
with numerical value approximately equal 0.4, u* G is the surface friction 
velocity, and j2f u (z/L 0 ) is a universal function of z/L 0 . The quantity L 0 
is the surface Monin-Obukhov stability length, namely. 


L 

o 


u fo C P P l 

k xS H o * 


(A-2) 


In this equation H Q is the surface heat flux, p and f denote the mean flow 
density and Kelvin temperature, g is the acceleration of gravity, and Cp 
is the specific heat at constant pressure. The dimensionless shear 0 U 
is related to the flux Richardson number through the experimentally 
derived relationship 


0 U - (1 - ISRi) 1 / 4 , 


(A- 3) 


which is given in reference 3. The flux Richardson number is defined as 


Ri 


s. de 

ft 


(du/dz) 2 


(A-4) 


where 0 is the mean potential temperature at height z. The flux Richardson 
number is a function of z/L 0 . We shall invoke the Businger hypothesis [8] 
to relate Ri to z/L 0 , so that 

Ri = z/L 0 . (A-5) 
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Upon combining equations (A-l), (A- 3) , and (A-5) and integrating the 
resulting relationship, we find that 


where 


u(z) 






(1 + 18^- x / 4 
I 


d|. 


(A- 6) 


(A- 7) 


We have used the condition that u(z 0 ) = 0 in the derivation of (A-6), 
where z G is the surface roughness length. Equation (A-7) can be evalu- 
ated numerically for any value of z/L q . The lower bound of this integral 
may be set equal to zero because the contribution to \|r from the region 
0 < -z/L 0 < -z 0 /L 0 is negligibly small. 

Equations (A-4) through (A-7) can be used to calculate the scaling 
velocity u* D and L 0 . The procedure for calculating these quantities is 
as follows: (1) calculate the gradient Richardson number, equation (A-4), 

with the mean flow wind and potential temperature profile data in the 
Monin layer, (2) calculate L Q with equation (A-5), and finally (3) cal- 
culate u* 0 with (A-6) and (A-7). 

The 18- and 30-meter temperature and wind data were used to estimate 
the Richardson number. The distributions of u(z) and 0(z) between the 
18- and 30-meter levels were assumed to be logarithmic profiles. The 
expressions for u(z) and 0(z) were differentiated with respect to z and 
evaluated at the geometric height z = 23 meters to yield estimates of 
bn/bz and bd/bz for the calculation of Ri. This Richardson number was 
used to calculate L 0 in step (2) above. The 18-meter level wind speed 
and i|r(18m/L 0 ) were used to calculate u^q in step (3). 

The surface roughness lengths z D that were associated with the NASA 
150-meter meteorological tower site and that were used in the calculation 
of u^ D are given in reference 4. 
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APPENDIX A 


Calculations of the Scaling Parameters u^ Q and L Q 


In the unstable Monin layer, the dimensionless mean flow shear is a 
universal function of z/L 0 , so that 


u *o d z 


0 u ( z/L o>’ 


(A-l) 


where u(z) is the mean wind speed at height z, k x is von Karman' s constant 
with numerical value approximately equal 0.4, u* 0 is the surface friction 
velocity, and j2f u (z/L 0 ) is a universal function of z/L 0 . The quantity L 0 
is the surface Monin-Obukhov stability length, namely, 


L 

o 


“lo V * 

k lS H o * 


(A-2) 


In this equation H 0 is the surface heat flux, p and T denote the mean flow 
density and Kelvin temperature, g is the acceleration of gravity, and Cp 
is the specific heat at constant pressure. The dimensionless shear 0 U 
is related to the flux Richardson number through the experimentally 
derived relationship 


0 U - (1 - 18Ri) l/4 , 


(A- 3) 


which is given in reference 3. The flux Richardson number is defined as 


Ri 


£L d9 

a & 

(du/dz) 2 


(A-4) 


where 0 is the mean potential temperature at height z. The flux Richardson 
number is a function of z/L 0 . We shall invoke the Businger hypothesis [8] 
to relate Ri to z/L q , so that 


Ri = z/L q . 


(A- 5) 
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Upon combining equations (A-l), (A-3), and (A-5) and integrating the 
resulting relationship, we find that 


u(z) 



where 


(A-6) 




a ± 18 I )- 3 -/ 4 

s 


d|. 


(A- 7) 


We have used the condition that u(z Q ) = 0 in the derivation of (A-6), 
where z G is the surface roughness length. Equation (A-7) can be evalu- 
ated numerically for any value of z/L Q . The lower bound of this integral 
may be set equal to zero because the contribution to from the region 
0 < -z/L 0 < -z 0 /L 0 is negligibly small. 

Equations (A-4) through (A-7) can be used to calculate the scaling 
velocity u* D and L 0 . The procedure for calculating these quantities is 
as follows: (1) calculate the gradient Richardson number, equation (A-4), 

with the mean flow wind and potential temperature profile data in the 
Monin layer, (2) calculate L Q with equation (A-5), and finally (3) cal- 
culate u* c with (A-6) and (A-7). 

The 18- and 30-meter temperature and wind data were used to estimate 
the Richardson number. The distributions of u(z) and 0(z) between the 
18- and 30-meter levels were assumed to be logarithmic profiles. The 
expressions for u(z) and 9(z) were differentiated with respect to z and 
evaluated at the geometric height z = 23 meters to yield estimates of 
du/dz and dd/dz for the calculation of Ri. This Richardson number was 
used to calculate L 0 in step (2) above. The 18-meter level wind speed 
and T|r(18m/L 0 ) were used to calculate u^ 0 in step (3). 

The surface roughness lengths z 0 that were associated with the NASA 
150-meter meteorological tower site and that were used in the calculation 
of u* 0 are given in reference 4. 
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APPENDIX B 


Wind Speed and Temperature Profile Data and Other Parameters 


The wind profile and temperature data that were used in this report 
are given in Tables B-l and B-2. The 18-meter level wind direction is 
tabulated here because the surface roughness length at the NASA 150-meter 
tower facility is a function of wind direction. The values of Ri(23 m), 
L q , and u* Q that were calculated with these data are given in Table B-3. 


TABLE B-l 

Table of Temperature Profile Data* 
(Temperature in °F) 


Case No. 

Date 

Time 

(EST) 

T(3m) 

(AT) x 

(AT) 2 

299 

1/23/68 

1315-1400 

76.0 

-1.41 

-1.73 

305 

1/26/68 

1130-1230 

52.4 

-2.19 

-2.66 

310 

2/8/68 

915-1015 

40.5 

- .32 

-1.82 

319 

2/26/68 

1030-1130 

59.0 

-1.28 

-2.63 

355 

3/20/68 

1330-1430 

75.0 

-1.82 

-2.21 

365 

3/28/68 

905-932 

73.5 

-1.44 

-2.42 

366 

3/28/68 

1145-1245 

74.5 

-2. 63 

-3.13 

406 

4/7/68 

1507-1537 

79.0 

-2.02 

-2.57 

445 

5/7/68 

1400-1415 

78.2 

-2.29 

-2.54 

551 

6/29/68 

947-1027 

79.8 

- .93 

-1.90 

554 

6/30/68 

933-952 

83 . 1 

-1.40 

-1.91 


*(AT) X = T(18m) - T (3m) . 
(AT) 2 = T (30) - T (3m) . 
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